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Abstract 

We show that the present limit from CLEO on the inclusive decay 6 ^ 57 pro- 
vides strong constraints on the parameters of the charged Higgs sector in two-Higgs- 
Doublet-Models. Only a slight improvement in the experimental bound will exclude 
the region in the Supersymmetric Higgs parameter space which is inaccessible to 
collider searches. 



In the Standard Model (SM) of electroweak interactions, spontaneous symmetry 
breaking is achieved through a single Higgs doublet, which then generates masses for both 
the gauge boson and fermion sectors. However, many extensions of the SM predict the 
existence of two Higgs doublets I]!]]. The physical spectrum of these models consists of three 
neutral Higgs scalars, two of which are CP-even {h^,H^) while one is CP-odd and 
two charged Higgs scalars (H^). In general, tree-level flavor changing neutral currents 
are present in two-Higgs-Doublet Models (2HDM), but can naturally be avoided if each 
fermion type receives mass from the vacuum expectation value (vev) of a single doublet. 
In one such model, hereafter referred to as Model I, one doublet (02) provides masses for 
all fermions and the other doublet (0i) decouples from the fermion sector. This model 
predicts distinctive phenomenology and has recently been revived in the literature0. In a 
second model (Model II), 02 gives mass to the up-type quarks, while the down-type quarks 
and charged leptons receive their masses from 0i. Supersymmetry and many axion theories 
predict couplings of the type present in Model II. 

Many different types of experiments may reveal the existence of physics beyond the 
SM. In the case of 2HDM, one could hunt for the enlarged Higgs spectrum directly in 
both e~^e~ and hadron colliders. The strongest direct search limits on these particles are 
provided by LEP^; in particular, the mass of the if^ is restricted to be > Mz/2 in all 
2HDM. It is also possible that such new physics may manifest itself indirectly in low-energy 
phenomena. Here, we will show that powerful bounds on the charged Higgs sector in 2HDM 
arise from the radiative b-quark decay, b 57, and that these constraints surpass those 
from direct collider searches. As we will see below, the application of these bounds to 
the particular case of supersymmetric models, could close the window of parameter space 
where there are no detectable Higgs signatures at colliders. A 90% C.L. upper bound on 
the branching fraction for this mode, B{b 57) < 8.4 x 10~^, has been obtained by the 
CLEO Collaboration via an examination of the inclusive photon spectrum in 5-meson 



1 



decays. Recent detector refinements coupled with increasing integrated luminosity leads us 
to anticipate that either the current limit will be strengthened, or the decay may actually 
be observed in the near future. 

Before presenting our results, we first briefiy describe the nomenclature of 2HDM. 
Each doublet obtains a vev fj, subject only to the constraint that v\ + = f ^, where v 
is the usual vev present in the SM. The charged Higgs interactions with fermions are then 
governed by the Lagrangian 

^ = o /^n^ H^[Vijmu,AuUi{l - -i^)dj (1) 

+ Vijmd^AdUi{l + 75)(ij + mi>Aii'{l + 75)^] + h.c. , 

where g is the usual 811(2)^ coupling constant and Vij is the appropriate element of the 
Kobayashi Maskawa (KM) matrix. In Model I, Au = cot (3 and A^ = Ae = — cot (3, and in 
Model II, Au = cot (3 and A^ = Ai = tan/?, where tan (3 = V2/V1 is the ratio of vevs. In 
a general 2HDM, the mass and tan/5 are a priori free parameters, as are the masses 
of all the neutral Higgs fields. However, in supersymmetric models, mass relationships 
exist between the various Higgs scalars. At tree-level, in such models, only two parameters 
are required to fix the masses and couplings of the entire scalar sector, but once radiative 
corrections are included 0, the values of the top-quark and squark masses also need to be 
specified. 

The transition b —>■ s'-f proceeds through electromagnetic penguin diagrams, which 
involve the top-quark, together with a or SM boson in the loop. The expression 
for the branching fraction in 2HDM is given in the literature To fully appreciate our 
results, it is important to observe the tan/5 dependence in the transition amplitude for both 
models. At the W scale the coefficients of the operators which mediate this transition take 
the generic form 

Ci{Mw) = Aw{ml/M^) + AA},(m?/m^±) + -^Ajj{m^Jmli) , (2) 
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where A = —1/ +1 in Models I and II, respectively. Aw corresponds to the SM 

12 -I- 

amplitude and Aj^ represent the H contributions; their analytic form is given for each 
contributing operator in Ref. 6. We employ the explicit form of the QCD corrections as 
stated in Eqns. (4.2) and (4.6) of Grinstein et al.^. In order to evaluate these corrections, 
we use the 3-loop expression for and fit the value of the QCD scale A to obtain consistency 
with measurements]^] of ^^(Ml) at LEP. In the SM, this procedure yields B{b —>■ 57) = 
(2.56 - 3.94) X 10"^ for rut in the range 90 - 200 GeV. 

The dependency of the branching fraction on m^, mfj^, and tan (3 in both Models I 
and II is explicitly presented in Ref. 7. Enhancements over the SM rate only occur in Model 
I for small values of tan j3. This is due to the fact that the contributions to the amplitude 
are always scaled as cot^ (3, as seen in Eqn. 2. We also note that a destructive interference 
between the and W contributions is possible for some values of the parameters due to 
the minus sign in Eqn. 2. In Model II, large enhancements also appear for small values of 
tan (3, but more importantly, the branching fraction is found to always be larger than that 
of the SM. This occurs independently of the value of tan/5 due to the presence of the term 
A]j. For certain ranges of the model parameters, the resulting value of 5(6 — S7) exceeds 
the CLEO bound, and consistency with this limit thus excludes part of the — tan jS 
plane for a fixed value of mj. This is shown in Fig. 1 for both models, where the excluded 
region lies to the left and beneath the curves. We see that in Model I only a small region 
of the parameter space is excluded, regardless of the value of rrit, however, in Model II 
a lower limit on m//± is obtained for all values of tan/? for the reasons discussed above. 
For example, with rrit = 150 GeV, we find that mH± > 110 GeV at large tan/?, and that 
even stronger bounds on mH± result for values of tan/? ^1. If the CLEO limit were to 
substantially improve, the constraints in Model II would strengthen drastically, while those 
in Model I would be essentially unchanged. In performing our calculations, we have tried 
to be as conservative as possible. We scale T{b ^ s-y) to the 6-quark semi-leptonic decay 
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width employing both QCD and phase space corrections (taking mc/mh = 0.3) to b ^ ceue, 
and take B{b cez/g) = 10.8% (Ref. 8). This procedure removes the sensitive dependence 
of the b sj rate to the values of the KM matrix elements. The precise value of rrify (we 
set rrLf, = 5 GeV) that is used in the QCD corrections and phase space factors can have a 
sizable effect. For example, if = 4.5 GeV were taken, the bound on mfj± in Model II 
would increase to 165 GeV for rrit = 150 GeV in the large tan/3 limit. QCD corrections 
which include the effects of separate mass scales for rrit > M\y are found to increase the 
radiative decay rate by a few percent, but have not been used here, since a separate scale 
for mH± should also be taken into account. However, we expect that this effect would not 
be too large since the running of as is slow between pairs of heavier mass scales. 

To show the power of this indirect search technique, we present in Fig. 2 the branch- 
ing fraction for Model II in the limit of large tan /5 as a function of with mH± = mt — mf,. 
If the actual value (or future upper bound) for the branching fraction were to lie below the 
solid curve, then the decay t — > bH^ is kinematically forbidden for a particular value of 
rrit. For example, the present CLEO result implies that t — > bH^ is already excluded for 
rrit ^ 215 GeV. We see that if the CLEO upper limit were to improve to 6 x 10~^, then 
this decay mode of the top-quark would be excluded for all values of the model parameters. 
Here, we made use of the large tan/3 limit, since it minimizes the contributions to 
b sj. 

In the supersymmetric case, the bounds shown in Fig. lb are more conventionally 
displayed as an allowed region in the tan/3 — rriA plane, where ruA is the mass of the 
CP-odd field. This is displayed in Fig. 3a for various values of rrit, where the radiative 
corrections to the SUSY mass relations have been employed assuming Msusy = 1 TeV. 
We find that our results depend only weakly on the exact value chosen for the squark 
masses. For rrit = 150 GeV, the excluded region is comparable to what can be explored IPH 



by LEP I and II. The region of supersymmetric parameter space which remains is exactly 
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that in which the hghtest CP-even scalar behaves hke the SM Higgs. Figure 3b shows 
the growth in the size of the excluded region if the CLEO bound were to improve to 
B{b S7) < 4, 5, 6, or 7 X 10-^ assuming = 150 GeV. If 5 < 6 x 10"^ then the 
window of parameter space left uncovered [p!0| ] by both e+e~ and hadron collider searches 
would be excluded. This would imply that the SSC/LHC could cover the entire remaining 
allowed supersymmetric Higgs parameter region. 

We note that in supersymmetric theories, other super-particles can also contribute 



to the one- loop decay 6 — > 57, and generally lead to a further enhancement in the rate[|ll 
However, a complete examination of the full supersymmetric parameter space needs to 
be performed to determine if supersymmetric contributions always enhance the branching 
fraction. If that were the case, then the limits presented in this paper could be strengthened 
in the supersymmetric version of the 2HDM. However, if some range of parameter values 
yielded a destructive interference in the 6 — > 57 amplitude, the bounds in the supersym- 
metric 2HDM could become weaker. 

In conclusion, we have shown that the decay 6 — > 37 is by far the most restrictive 
process in constraining the parameters of the charged Higgs sector in 2HDM, yielding 
bounds which are stronger than those from other low-energy processes and from direct 
collider searches. We anxiously await future results from the CLEO experiment. 
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Figure Captions 



Figure 1. The excluded regions in the m//± — tan/? plane for various values of rut, resulting 
from the present CLEO bound in (a) Model I and (b) Model II. In each case, from 
top to bottom, the sohd (dashed dot, sohd, dotted, and dashed) curve corresponds to 
rrit = 210(180, 150, 120, and 90) GeV. The excluded region lies to the left and below 
each curve. 

Figure 2. B{b — > S7) as a function of rrit, with m^± = mt — rrih in the large tan (3 limit in Model 
II. 

Figure 3. (a) The excluded region from the present CLEO limit in the tan — plane for 
various values of rrit ^ indicated, (b) Excluded regions for rrit — 150 GeV, if the 
CLEO bound was improved to i? < 4, 5, 6 or 7 x 10"*^. In each case, the excluded 
region lies to the left and below each curve. 
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